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Abstract 

We study trace functions on the form t — > TV f(A + tB) where 
/ is a real function defined on the positive half-line, and A and B 
are matrices such that A is positive definite and B is positive semi- 
definite. If / is non-negative and operator monotone decreasing, then 
such a trace function can be written as the Laplace transform of a 
positive measure. The question is related to the Bessis-Moussa-Villani 
conjecture. 

Key words: Trace functions, BMV-conjecture. 

1 Introduction 

The Gibbs density matrix of a system with Hamiltonian H in equilibrium 
and temperature kT = 1//3 is given by exp(— (3H). One would like to study 
perturbations H + XHi of an exactly solvable Hamiltonian H and see how 
the thermodynamical quantities are changed. This question was studied by 
Bessis, Moussa and Villani in [2] where it is noted that the Pade approximant 
to the partition function Z(f3) = Tr exp(— (3(H + XHi)) may be efficiently 
calculated, if the function 

A^Trexp(-/3(iJo + A#i)) 

is the Laplace transform of a positive measure. The authors then noted that 
this is indeed true for a system of spinless particles with local interactions 
bounded from below. The statement also holds if H and Hi are commuting 
operators, or if they are just 2x2 matrices. These observations led to the 
formulation of the following conjecture: 
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Conjecture (BMV). Let A and B be n x n matrices for some natural 
number n, and suppose that A is self-adjoint and B is positive semi- definite. 
Then there is a positive measure \i with support in the closed positive half-axis 
such that 



for every t > 0. 

The Bessis-Moussa-Villani (BMV) conjecture may be reformulated as an 
infinite series of inequalities. 

Theorem (Bernstein). Let f be a real C°° -function defined on the positive 
half-axis. If f is completely monotone, that is 



then there exists a positive measure ji on the positive half-axis such that 



for every t > 0. 

The BMV-conjecture is thus equivalent to saying that the function 



is completely monotone. A proof of Bernstein's theorem can be found in [3]. 

Assuming the BMV-conjecture one may derive a similar statement for 
free semicircularly distributed elements in a type Hi von Neumann algebra 
with a faithful trace. This consequence of the conjecture has been proved by 
Fannes and Petz [6J. A hypergeometric approach by Drmota, Schachermayer 
and Teichmann jlj gives a proof of the BMV-conjecture for some types of 
3x3 matrices. 

1.1 Equivalent formulations 

The BMV-conjecture can be stated in several equivalent forms. 
Theorem 1.1. The following conditions are equivalent: 

(i). For arbitrary n x n matrices A and B such that A is self-adjoint and 
B is positive semi-definite the function f(t) = Trexp(A — tB), defined 
on the positive half-axis, is the Laplace transform of a positive measure 
supported in [0, oo). 




( _ ir/ (n) (t) > t>0,n = 0,1,2,..., 




fit) = Tr exp(A-tB) 



t > 
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(%%). For arbitrary n x n matrices A and B such that A is self-adjoint and 
B is positive semi-definite the function g(t) = Tr exp(A + itB), defined 
on the positive half-axis, is of positive type. 

(Hi). For arbitrary positive definite n x n matrices A and B the polynomial 
P(t) = Tr(A + tB) p has non-negative coefficients for any p = 1, 2, . . . . 

(iv). For arbitrary positive definite n x n matrices A and B the function 
ip(t) = Tr exp(A + tB) is m-positive on some open interval of the form 
(—a, a). 

The first statement is the BMV-conjecture, and it readily implies the 
second statement by analytic continuation. The sufficiency of the second 
statement is essentially Bochner's theorem. The implication (Hi) =^> (i) is 
obtained by applying Bernstein's theorem and approximation of the expo- 
nential function by its Taylor expansion. The implication (i) =>- (Hi) was 
proved by Lieb and Seiringer [TH]. A function ip : (—a, a) — ► R is said to be 
m-positive, if for arbitrary self-adjoint k x k matrices X with non-negative 
entries and spectra contained in (—a, a) the matrix <p(X) has non-negative 
entries. The implication (Hi) ^> (iv) follows by approximation, while the 
implication (iv) =3- (i) follows by Bernstein's theorem and [HI Theorem 3.3] 
which states that an m-positive function is real analytic with non-negative 
derivatives in zero. 



2 Preliminaries and main result 

Let / be a real function of one variable defined on a real interval /. We 
consider for each natural number n the associated matrix function x — > f(x) 
defined on the set of self-adjoint matrices of order n with spectra in /. The 
matrix function is defined by setting 

v v 

f(x) = Y,f( Xi ) Pi where x = J2 X * P * 

i=l i=l 

is the spectral resolution of x. The matrix function x — > f(x) is Frechet 
differentiable [2] if J is open and / is continuously differentiable [3]. The 
norm of the Frechet differential df(x) may be an unbounded function of the 
order m, cf. [1^1 [THJ . If however / is assumed to be twice differentiable, then 
the norm of df(x) is bounded independently of the order n for all self-adjoint 
matrices x with spectra contained in a fixed compact subset of J, cf. 
Corollary 2.9] and [12] • We consider in this article the function 

ip(x) = Trf(x) 
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defined on the set of self-adjoint matrices of order n with spectra in /. The 
Frechet differential is given by dip(x) = Tr df(x), cf. [10J. 



2.1 The BMV-property 



Definition 2.1. A function /: R + — > R is said io /«ai>e ifte BMV-property, 
if to each n = 1, 2, . . . and eae/i pa«r of nx n matrices A and B, such that A 
is positive definite and B is positive semi-definite, there is a positive measure 
H with support in [0, oo) such that 



for every t > 0. 

The BMV-conjecture is thus equivalent to the statement that the function 
t — ► exp(— t) has the BMV-property. The main contribution of this paper is 
the following result. 

Main Theorem. Every non-negative operator monotone decreasing function 
defined on the open positive half-line has the BMV-property. 

3 Differential analysis 

An simple proof of the following result can be found in [TJl Proposition 1.3]. 

Proposition 3.1. The Frechet differential of the exponential operator func- 
tion x — > exp(x) is given by 



where A(s) = exp(sx)hexp(—sx) for s 6 R. 

This is only a small part of the Dyson formula which contains formalisme 
developed earlier by Tomonaga, Schwinger and Feynman. The subject was 
given a rigorous mathematical treatment by Araki in terms of expansionals 
in Banach algebras. In particular [TJ Theorem 3], the expansional 





exp(sx)h exp((l — s)x) ds 




E r (h;x) = / / •••/ A(s n )A(s n . 1 )---A(s 1 )ds n ds. 



n-1 



■ ■ ■ ds 



i 



Jo Jo Jo 



is absolutely convergent in the norm topology with limit 



E r (h; x) = exp(x + h) exp(— x). 
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We therefore obtain the pth Frechet differential of the exponential operator 
function by the expression 



d p exp(x)h 



V 



1 PS 



>p— 1 



= p\ / / ••• / A(sp)A(s p -i) ■ ■ ■ A(si) exp(x) ds p dsp-i- ■ -dsi. 
Jo Jo Jo 

3.1 Divided differences 

The following representation of divided differences is due to Hermite [I 
confer also [TTj fill]. 

Proposition 3.2. Divided differences can be written in the following form 
[x ,xi]/ = J f'((l - t^xo + *ia:i) dt 
[x ,x 1 ,x 2 )f = J J f"(\A-ti)xo + {t 1 -t 2 )x 1 + t2X^jdt2dti 



[x ,x ir -- ,x n ] f =1 [ - I f {n) ({l-t 1 )x + {t 1 -t 2 )x 1 + 



JO 



+ (fn-1 ~ tn) x n-l + ^n^n^ dt n ■ ■ ■ dt 2 dt\ 



where f is an n-times continuously differential function defined on an open 
interval I, and xq, xi, . . . , x n are (not necessarily distinct) points in I. 

Lemma 3.3. For k = 1, 2, . . . and real numbers Xi, . . . , A& we have 

[ e~ MS [sAi, . . . , sX k ] cxp s^ 1 ds = t k e~^[tX h ...,t\ k , tfi} exp 
Jo 

for any real [l and t > 0. 
Proof. For k = 1 we calculate 

I e-^{s\ l ] cxp ds= [ e s ^-tids = te- fit [tXi,ti4exp- 
Jo Jo 
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Assuming the formula valid for k we obtain for k + 1 the expression 

/ e-^[sX 1 ,...,sX k+ i\e X pS k ds 
Jo 

f l c -^ s [ S -^li • • • ; S ^fc]cxp ~ [sX 2 , ■ ■ ■ , sA fc+ l]cxp s k ^ 

Jo sXi — sXk+i 

= t ~T (^e-^tXu tX k , t/i] exp - t k e-^[tX 2 , tX k+1 , tfi] cxp ) 



\fe+l 

t fc+1 e-^[tA 1 ,...,tA fc+1 ,t/i] 



cxp 



provided Ai 7^ A^+i. The case Ai = A^+i then follows by continuity, and the 
lemma is proved by induction. QED 

Theorem 3.4. Letx andh be operators on a Hilbert space of finite dimension 
n written on the form 

n n 

x = Xita and h = h^e^ 

i=l i,j=l 

where {ejj}™- =1 is a system of matrix units, and Ai, . . . , X n and hij for i,j = 
1, . . . , n are complex numbers. Then the pth derivative 

dP nn / , X I 

-— Tr expix + th)\ 
dtv FV J \t=o 

n n 

^ ^ ' ' ' ^ ^ tl'ipip-i ' ' ' ^1211 hi\ip [^H ' 1 ' ' ' 1 ^i P i ^jplcxpj 
il=l i P =l 

where [A^, Aj 2 , • • • , Aj p , AjJ exp are divided differences of order p + 1 of the 
exponential function. 

Proof. We first prove the formulas 

4(so) = / / •••/ A(s k )A(s k -i) ■ ■ ■ A(s 1 )ds k ds k - 1 - ■ -ds! 
Jo Jo Jo 

. . n n n 

I 1 ) s fc W • • • V /; /, h- h 

3=1 h=l «fc=l 

X e S()Aj [So-^in • • • , S()Xi k , S Xj} cxp £'i k j 

for natural numbers k = 1, . . . ,p and real numbers sq > 0. For k = 1 we 
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calculate the integral 

h(so) = / A(si) dsi = I exp(six)hexp(—six)dsi 



o 



^ ^ / exp(six)e iu exp(-six) <isi 

i 1= l i=l J ° 

Y Y / e M s i(K - \))e n j 

U=l 7=1 J ° 



11=1 J 

n n 



exp Giij 

3=1 li=l 

in accordance with (JXJ) . For k > 2 and assuming the formulas ((H) valid for 
k — 1 we obtain the expression 

■ ■ ■ / A(s k )A(s k - 1 ) ■ ■ ■ A(s 1 ) ds k --- ds 2 ds 1 
Jo 

so 



J fc _i(si)A(si) ds 



i 



rs n 

/ V ft ft- ■ ft 

/ °1 / ■ l h k i k _ 1 "■1312"' 

JO . . , 



x [siAj 2 , . . . , Si\i k , siA m ] exp ei km A(s\) ds\. 



We then insert 



= Y Y ex P( s i( A u - A i)) e ni 

and using e ikm e ixJ - = 5(m,ii)e ik j and the symmetry of the divided difference 
we obtain the expression 



^ hi k i k _ 1 ' ' ' h>i 2 ii I 6 J [siAjj , SiAj 2 ) • • • ; >§i Aj fe ] ex p S;L ds% 6j fe 

J,»i,i2,...,tfc=l 

for Ik(so). Finally, using Lemma EP1 we calculate 

n 

Ik(so) = $0 ^ ^ hi k i k _ 1 ■ ■ ■ h{ 2 i 1 Hi^jC 3 [soAi i; . . . , SoAj fe , SoAj]exp CiiU 
j,ii,ia,...,tfc=l 

which proves ((H) by induction. We next observe that 

d p exp(x)h p = p\ I p (l) exp(x) 
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where differentiation is with respect to x when nothing else is indicated. 
Finally, since 

dP 



dtP 

we obtain 



exp(x + th) = d p x+th exp(x + th)h p 



■ Tr exp(a; + th) | = Tr (d p exp(x)h p J = p\ Tr (j p (i) exp(a;^ 



p\ Tr 



^ ^ hipip-i ' ' ' hi 2 i 1 hi 1 je J [Ajj, . . . , Aj p , Aj] exp &% v j exp(x) 



j,il,i2,...,i p =l 
n 



which is the statement of the theorem. QED 

Lemma 3.5. Let f G C 00 (R), and let x and h be self-adjoint operators on 
a (possibly infinite dimensional) Hilbert space H. Then the operator function 
x — > f(x) is infinitely Frechet dijferentiable and the pth Frechet differential 
is for p > 1 given by 



d p f{x)h p = / d p x exp(-isx)h p f(s)ds, 

J — oo 

where d x indicates differentiation with respect to x and 



i r°° 

f(s) = — / exp(tts)f(t)dt 



is £/ie Fourier transform of f. 

Proof. We note that the statement is true for p = 1 by (T2J Theorem 1.5] 
and assume the statement of the lemma to be valid for p. It follows from the 
definition of the Frechet differential that the expression 

d p x+h exp(— is(x + h))h p — d p exp(—isx)h p — d x (d p exp(—isx)h p )h 

even after division by \\h\\ tend to zero as h — > 0. We then multiply the above 
expression by the Fourier transform / and integrate. By Lebesgues's theorem 
of dominated convergence we therefore obtain that also the expression 

/oo 
d x (d p x exp(-isx)h p )hf(s) ds, 
-oo 
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even after division by \\h\\, tend to zero as h — > 0. Hence 

/oo 
< +1 exp(-isa;)/i p+1 /(s) ds 
-oo 

and the lemma is proved by induction. QED 
In the next corollary we need the identity, 

(2) t p - 1 [tX 1 ,...,t\ p ] f = [X 1 ,...,X p ] ft where f t (s) = f(ts), 

valid for p times continuously differentiable functions /. The statement is 
easily proved by induction after p. 

Corollary 3.6. Let f : I — > R be a C°°-function defined on an open and 
bounded interval I, and let x and h be self-adjoint operators on a Hilbert 
space of finite dimension n written on the form 



n 



= XiCu and h = hi 
i=i «,i=i 



where {eij}fj =1 is a system of matrix units, and Ai, . . . , X n are the eigenvalues 
of x counted with multiplicity. If the spectrum of x is in I, then the trace 
function t — * Tr f(x + th) is infinitely differentiable in a neighborhood of zero 
and the pth derivative 

d p l 
— Trf(x + th)\ 
dtP Jy J \t=o 

n n 

= ^ ^ ' ' ' ^ ^ hi^hi^ ■ ■ ■ hi p l i p hi p i 1 [Ajj , Aj 2 , , Aj p , Ajjj, 

ii=l i p =l 

where [A^, Aj 2 , • • • , Aj p , AjJ/ are divided differences of order p+ 1 o/ £ne /unc- 
^on /. 

Proof. Since the spectrum Sp(x) is compact, there is an open and bounded 
interval J such that 

Sp(x) cJcJcI 
and we may extend the restriction f\j to a function g G C 00 (R). Since the 
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spectrum Sp(x + th) is contained in J for small t we obtain 



dP 

— Trf(x + th) 
dtP Jy ' 



d p 

t=o dtP 



t=o 



f°° dP 

I —— Tr exp(—isx — isth)) g(s)ds 
J-oo dt p t=o 

/qq n n 

x [-isX h ,--- ,-is\i p ,-is\ ip } exp g(s)ds 

n n poo 

• ^ ^ ' ' ' ^ ^ hi p i p _ 1 ■ ■ ■ hi^i x hi x i p I , \i p , Ai p ]exp(— is-) di s ) ds 

i 1= l i p =l 



ii=l «r,=l 



n n 



£*• ^ ^ ■ ■ ■ ^ ^ hi p i p _ 1 ■ ■ ■ hi^hi^ [At 1( , Aj p , Ajp]^ c?s, 

*x=l ip=l 

where g is the Fourier transform of g, and we used ([2|) and the linearity in g 
of arbitrary divided differences [Ai, . . . , \ p ] g . The statement now follows by 
noting that / = g in a neighborhood of the spectrum of x. QED 

4 Proof of the main theorem 

Proposition 4.1. Consider for a constant c > t/ie function 

1 



c + t 



t > 0. 



For arbitrary n x n matrices x and h such that x is positive definite and h 
is positive semi-definite we have 



dP 

:-iy—Trg(x + th) 



> 



t=o 



for p= 1,2, .... 

Proof. We first note that the divided differences of g are of the form 

(3) [A l5 A 2 , . . . , X p } g = (-ir- 1 g(X 1 )g(X 2 ) ■ ■ ■ g(X p ) p = 1, 2, . . . 

for arbitrary positive numbers Ai, A 2 , . . . , X p . There is nothing to prove for 
p = 1. Assume the statement true for p > 2 and notice that 

g{\)-g(») (c + A)- 1 - (c + zi)- 1 



X — ji 



X — n 
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for A 7^ /i. Therefore 

[Ai, . . . , A p , A p+ i] 



[Ai, . . . , \p]g — [A 2 , . . . , A p+ i] g 
Ai — A p+1 



= (-ir L 9(x 2 )---9(x P y- x , 

= (-l)^(Ai).-^(ApMVi) 

for Ai 7^ A p+ i and the general case follows by approximation. Next, applying 
Corollary 13.61 we obtain 

d p I 
— lTg(x + th)\ 
at p I t=o 



n n 



«i=l ip=l 

n n 



= (-i) p p\ J2"'Y1 h ™Ki 3 ■ ■ ■ K-iipKn g(K)g(K 2 ) • • -g(K)9(K)- 

ii=l i p =l 

Since h is positive definite, it is of the form h = aa* for some matrix a and 
therefore 

n n 

hij ^ ^ O'im^'mj ^ ^ Q'imQ'jm (^i | ^; J 1, . . . , 71, 

m=l m=l 

where 

di = ( Oil • ■ ■ a in ) G C" 

is the ith row in the matrix a. We then set £j = g(Xi)ai and 6j = g(Aj) 1//2 cij 
for i = 1, . . . , n and calculate 



— Trg(x + th) 



p\ dt p 



n n 



t=0 



^ ' ' ' E ^1*2^2*3 ' ' ' hip-iiphiph 9(^1)9(^12) ' ' '#(^ip)#(<^il) 

«i=l ip=l 
n n 

Yl ' ' I h n)(K 1 6»s ) • • • (K-i I & i P )( & i P I 



il=l i p =l 

But any sum of the form 



n n 



12=1 i P =l 
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is non-negative. It is obvious for p = 2 since then 



»2=1 



and for p = 3 since then 



EE(£ I 6 «)( 6 « I 6 «)( 6 « 10= E^ I 6«)6i 2 I $^(0 6 is )6i 

»2=1 »3=1 V2 = l »3=1 

For p > 4 we may therefore use induction over either the even or the odd 
natural numbers by noting that 

n n 

s = E • • • i i • • • (Vi i i o 

12=1 i p =l 
n n 

= E " ' E fa I b is)( b is \K)--- (K-2 1 I ^) 

13=1 ?p-l = l 

where the vector 

n 

v = E(£ i 

i=l 

This concludes the proof of the statement. QED 
Proof of the main theorem. Consider again the function 

for c > and arbitrary nxn matrices x and h such that x is positive definite 
and h is positive semi-definite. We first note that 



dtp 



d p 

Trg(x + th) = — Tr g(x + t h + eh) 



t=to 



deP 



e=0 



for p = 1,2,... and to > 0. The function t — > Trg(x + th) is therefore 
completely monotone. Let now /: R + — > R be a non-negative operator 
monotone decreasing function. Any operator monotone decreasing function 
defined on the open positive half-line, thus in particular the function /, is 
necessarily of the form 



(4) 



f(t) = at + p + f (^- t ~) MO, 
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where a < and v is some positive Borel measure with support in [0, oo) 
such that the integral J(c 2 + l) -1 du(c) is finite, cf. |U Chap. II Theorem 1 
and Lemma 2]. Note that we may write 

r°° 1 — r f 
du(c) = / -■(c 2 + l)- 1 du(c) 



c + t c 2 + 1 / J c + t 

where for each t > the function c — > (1 — ct)(c + t)^ 1 is decreasing and 
bounded between —t and t" 1 . Since / has a finite limit for t — > oo one may 
derive that a = 0, and by appealing to Lebesgues's theorem of monotone 
convergence we also obtain f c{c 2 + l)" 1 dv{c) < oo, hence / allows the 
representation 

f(t)=/3+ rj—du( C ), 



jo c + t 

where we have incorporated the constant contribution from the integral in 
(|1J) into P such that 

P = lim /(t) > 0. 

t^oo 

We conclude that the function t — > Tr /(x + t/i) is completely monotone and 
thus by Bernstein's theorem is the Laplace transform of a positive measure 
with support in [0, oo). QED 

4.1 Further analysis 

One may try to use the Hermite expression in Proposition 13.21 to obtain a 
proof of the BMV-conjecture. Applying Theorem 13.41 and calculating the 
third derivative of the trace function we obtain 



Id 3 n 

- — Trexp(x - th)\ _ = (a p \ a^fai | | a p )[X p XiXjX p ) exp 

it— o . . 



3! dt 3 



f f f 

= / / / ( a p\ a i)i a i I a>j){aj | Op)exp((l - (ti -t 3 ))X p 

Jo Jo Jo P)iJ=1 

+ (ti - t 2 )Xi + (t 2 - t 3 )Xj) dt 3 dt 2 dt x 

where h = aa* and a* is the ith row in a. Assuming the BMV-conjecture this 
integral should be non-negative, and this would obviously be the case if the 
integrand is a non-negative function. 

Example 4.2. If we evaluate the above integrand in t\ = t 2 = 1 and t 3 = 1/3 
we obtain 



S = (a p | a,i)(ai | cij)(aj | a p ) exp ^ p — ^ dt 3 dt 2 dt\ 



a p | ajJ(Oj | a,J(aj | <2pjexp(^ 

p,»>j=i 
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If we in addition set 



ai =(1000,-10,1) 

a 2 = (-10, 10000, 1000) 

a 3 =(1,1000,202139) 

and (Ai, A2, A3) = 3 log 2 ■ (23, 11, 0), then the sum is an integer with value 

S = -487062506352658941731358505750. 

The values of S are extremely sensitive to the chosen figures, and they tend 
to be overwhelmingly positive. If for example the third entry in 03 is changed 
from 202139 to 202138, then 

S = 376189230591238013538921396773. 

The result is equally sensitive to changes in the values of (Ai, A 2 , A3). 

The above example has not been found by simulation. In fact, millions 
of simulations with randomly chosen entries have been carried out, without 
ever hitting a negative value of S. Instead, the example has been obtained 
by the study and proper modification of an example in [14J of two positive 
definite 3x3 matrices A and B such that Tr(BABAAB) = -3164. 

Another way forward would be to examine the value of loops of the form 

(ai I a 2 )(a 2 | a 3 ) • • • (a p _i | a p )(a p \ ax) 

since they, apart from an alternating sign, are the only possible negative 
factors in the expression of the derivatives of the trace functions. The value of 
a loop is a homogeneous function of degree two in the norm of the vectors, so 
we only need to consider unit vectors. The value of such a loop is furthermore 
invariant under unitary transformations and thus takes the minimal value in 
any subspace of dimension p. By applying a variational principle the lower 
bound 

- cos n ^— J < (ax I a 2 )(a 2 | a 3 ) • • • (a p _i | a p )(a p \ ax) 

was established in [TT|. The lower bound converges very slowly to —1 as p 
tends to infinity, and it is attained essentially only when all the vectors form 
a "/an" in a two-dimensional subspace. 

Remark 4.3. // we only consider one- dimensional perturbations, that is if 
h = cP for a constant c > and a one-dimensional projection P, then h is 
of the form h = {£i£j)i,j=i,...,n f or some vector £ = (£1, . . . , £ n ) and each loop 

h% Y i 2 hi 2 i z ■ ■ ■ hi p l i p hi p i 1 = H^iJI ■ ■ • ||£jp|| 



14 



is manifestly real and non-negative. This implies that the trace function 

t -> Trex.p(-(x + th)), 

for any self-adjoint n x n matrix x, is the Laplace transform of a positive 
measure with support in [0, oo). The same statement holds, with x positive 
definite, for the trace function t — > Tr f(x + th) associated with an arbitrary 
completely monotone function f. 
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